REPORT 1279 


THEORETICAL ANALYSIS OF INCOMPRESSIBLE FLOW THROUGH A RADIAL-INLET 
CENTRIFUGAL IMPELLER AT VARIOUS WEIGHT FLOWS 1 

By James J. Kramer, Vastly D. Pbian, and Chttng-Htta Wtj 


SUMMARY 

A method for the solution of the incompressible nonviscous 
ilow through a centrifugal impeller, including the inlet region, 
is presented. Several numerical, solutions are obtained for four 
weight flows through an impeller at one operating speed. These 
solutions are refined in the leading-edge region. The results are 
presented in a series of figures showing streamlines and relative 
velocity contours. A comparison is made with the results ob- 
tained by using a rapid approximate method of analysis. 

INTRODUCTION 

la order to provide the fundamental information about 
the internal flow necessary for the rational design of efficient 
centrifugal compressors, two-dimensional solutions of the 
potential flow through centrifugal compressors on both blade- 
to-blade and meridional surfaces are obtained by means of 
relaxation methods in references 1 to 4. In addition, a 
three-dimensional potential-flow solution is obtained by sim- 
ilar means in reference 5. However, all these solutions are 
for impellers with inducer sections extended infinitely far 
upstream or to the axis of the impeller and thus yield no 
information concerning the flow behavior ahead of and at 
the entrance to blades of finite thickness or blades which 
are not alined with the inlet stream. 

A rapid approximate method is developed in reference 6 
that will predict the blade surface velocities in centrifugal 
impellers. The accuracy of this method is investigated 
therein for the region downstream of the inducer section. 
However, because of the lack of an exact solution of the 
flow in the inlet region, there has been no verification of the 
approximate method in that region. 

Consequently, a method for analyzing the flow in a cen- 
trifugal impeller, including the inlet region, by numerical 
solution of the partial differential equation governing the 
flow was developed at the NACA Lewis laboratory and is 
presented herein. This method was applied to a 48-inch- 
diameter centrifugal impeller for both design and off-design 
flow conditions. The impeller is similar to that discussed 
in references 7 and 8. The results of the exact solution are 
compared with the results of the rapid approximate method 
of reference 6. 


In addition, the flow in the region of the leading edge of 
the blade was determined in more detail than that obtained 
in the original solution. This refinement of the original 
solution is obtained by solving the partial differential equa- 
tion in that region by relaxation methods using a grid of 
much finer mesh. More detailed information concerning the 
flow behavior near the leading edge than that obtained in 
the original solution is desirable, because a knowledge of 
the velocity gradients in this region is helpful in avoiding 
boundary-layer separation caused by rapidly decelerating 
flow. The magnitude of losses caused by boundary-layer 
separation at the leading edge is discussed in references 9 
and 10 for sharp-nosed blades, but the manner in which the 
losses occur is not discussed. 

ANALYSIS 

The formulation of the problem and the proposed method 
of solution of the problem are discussed in this section. 

STATEMENT OF PROBLEM 

The basic assumptions which are made concerning the 
physical nature of the flow determine the partial differential 
equation governing the flow. The assigning of proper 
boundary conditions to the problem then determines the par- 
ticular solution of the partial differential equation. 

Assumptions. — The flow is assumed to be steady, in- 
compressible, and nonviscous. The assumption of steady 
nonviscous flow is customary in compressor flow analyses. 
Several solutions have been obtained taking compressibility 
into account (e. g., refs. 1 and 2) . The addition of the con- 
dition of compressibility complicates the solution by a factor 
which was considered to be out of proportion with its value 
in this case. Therefore, the fluid was assumed incom- 
pressible. 

The further assumption is made that the flow is constrained 
to a blade-to-blade surface of revolution which is symmetrical 
about the impeller axis. Although the flow is constrained to 
this surface, a variation in the thickness of the stream sheet 
provides a closer approximation to the actual case. The 
shape of the stream surface in the axial-radial plane as well 
as the thickness variation are defined as functions of radial 
position, which are specified at the beginning of the solution. 


i Supersedes NAOA TN 3+48, ‘'Theoretical Analysis of Incompressible Flow Through a Radial-Inlet Centrilugal Impeller at Various Weight Flows. I— Solution by a Matrix Method 
and Comparison with on Approximate Method,” by Vasily D. Brian, James J. Kramer, and Ohung-Hua Wu, 1955; and TN 3«9, “Theoretical Analysis of Incompressible Flow Through a 
Radial-Inlet Centrifugal Impeller at Various Weight Flows, n— Solution In Leading-Edge Region,” by James J. Kramer. 1955. 
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The rear stagnation point is assumed to be located at the 
blade tip. The Kutta condition states that for a noncuspi- 
date blade with a sharp trailing edge the rear stagnation 
point occurs at the tip. However, for an impeller with a 
rounded trailing edge the location of the rear stagnation 
point cannot be predicted. It is necessary, therefore, to 
assume the location of the rear stagnation point. 

For the refined solution it is assumed that the values of the 
stream function obtained in the original solution along the 
boundaries of the region in which the refined solution is ob- 
tained remain unchanged during the numerical procedure of 
obtaining the refined solution. 

Differential equation. — In this analysis, the cylindrical 
coordinates r, 0, and z (see figs. 1 and 2) are used. (All sym- 
bols are defined in the appendix.) The angular velocity of 
the impeller is denoted by <o and the fluid density by p. The 
stream-sheet thickness in the 2 -direction is represented by b. 
The trace of the stream surface in the axial-radial plane is 
given by specifying 2 as a function of r. The slope drjdz of 
this curve is equal to the tangent of the angle between the 
axis of rotation and the tangent to the trace of the stream 
surface in the axial-radial plane and is denoted by X (see 
fig. 2) . Thus, the resultant velocity w is given by 



Axis of rotation 


Fiatnaz 2. — Axial-radial plane view of 48-inch-diameter centrifugal 

impeller. 

vd‘=W i e +W a r (l+^ 

where the subscripts r and 0 indicate components in the r- 
and ^-directions, respectively. 

The stream function Sf' is defined by the following differ- 
ential equations: 

-JL=_&pTi, e (la) 

(lb) 

In this report all derivatives with respect to r shall be under- 
stood to mean derivatives with respect to r on the stream sur- 
face; that is, 5 [dr in this report shall correspond to the bold- 
faced S/dr of reference 1 1, in which the differential equation 
for the type flow considered herein is derived. With those 



ANALYSIS OF INCOMPRESSIBLE FLOW THROUGH A CENTRIFUGAL IMPELLER AT VARIOUS WEIGHT FLOWS 651 


definitions and assumptions, the differential equation of the 
flow becomes (see ref. 11, p. 35) 

This equation, together with the boundary conditions, math- 
ematically determines the problem. 

Boundary conditions. — This analysis of the flow is a 
boundary-value problem of the first kind or a Dirichlet prob- 
lem. Certain boundaries of the flow and the values of the 
stream function on these boundaries are specified. Further- 
more, the flow is assumed to vary periodically in the circum- 
ferential direction, completing a cycle in one pitch angle, the 
angular distance between two adjacent blade mean lines. 
The rotational speed of the impeller and the weight flow 
through the compressor are alsp specified. The domain of 
the solution is extended sufficiently far upstream and down- 
stream so that the flow is assumed to be uniform at the up- 
stream and downstream boundaries. With the addition of 
these conditions, the problem is determined mathematically. 

For the refined solution the boundary conditions are ob- 
tained from the original solution. The value of the stream 
function, expressed as a dimensionless ratio of the weight 
flow through a single passage, on the blade surface is equal 
to 1 ; the values along the other boundaries are obtained from 
cross plots of the original data. The error involved in read- 
ing values from the cross plots was less than 0.0003 with 
values of 'i'/M ranging from 0 to 1. 

METHOD OF OBTAINING SOLUTION FOB FLOW THROUGH ENTIRE BLADE 

PASSAGE 

Superposition of four basic solutions. — The differential 
equation (eq. (2)) was solved by a superposition of four basic 
solutions. These four basic solutions form a set of linearly 
independent solutions such that all possible flows (including 
all tip speeds) are expressible as linear combinations of these 
basic solutions. The first of these, designated is a solu- 
tion of equation (2) with the condition that no flow crosses 
the upstream and downstream boundaries and co=wo?^0. 

The other three basic solutions, designated fa, fa, and fa, 
are solutions of the linear homogeneous equation obtained 
by equating the left side of equation (2) to zero. Thus, if 

r _ b! , /l - Aq-IA &I 

dr* r \r br /dr'VV "Tvse* 
then 4 * 1 , fa, and fa are solutions of 

L(fa= 0 

Because L is a linear operator, fa plus linear combinations of 
i/'i, fa, and fa will satisfy equation (2) for <o=«o. 


Boundary conditions for four basks solutions. — The flow 
region is represented by abcdefgh in sketch (a): 


e 



The upstream and downstream boundaries, ah and de, 
respectively, are placed sufficiently far from the blades so 
that flow conditions can be assumed uniform at these stations. 
The angular distance from a to h and from d to e is one pitch 
angle. For all the basic solutions, the condition that the 
flow is periodic about the axis of rotation with a period of 
one pitch angle makes it possible to obtain the solutions 
without a knowledge of the stream function along ab and gh. 
The finite-difference equation for points along these lines is 
obtained in the same manner as in reference 12. For the 
solution fa in which the flow is that induced only by the 
rotation of the impeller without any through flow, the value 
of V' along ah and de is specified zero, indicating no flow 
crossing the upstream and downstream boundaries. The 
values along the blade surfaces be and gf are also specified 
as zero. The solution to equation (2) for these boundary 
conditions is designated fa. 

The through-flow solution, that is, the flow through the 
stationary blade row, is obtained from linear combinations of 
three basic solutions to the homogeneous equation obtained by 
equating the left side of equation (2) to zero. All possible 
flows through the stationary blade row can be represented 
by linear combinations of these basic solutions. It can be 
seen from the boundary conditions shown in the following 
table that these three basic solutions are linearly independent: 


Basic 

solution 

Boundary values at point — 

a 

d 

e 

h 

to 

0 

0 

0 

0 

fa 

0 

0 

1 

1 

tl 

-1 

0 

1 

0 

fa 

—1 

-1 

0 

0 


The value of ^ along the trailing face is specified as zero and 
along the driving face as 1 for all three solutions. That 
these independent solutions are sufficient for the construction 
of all possible through flows can be seen from the following 
consideration. 
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Any flow is determined when the upstream and down- 
stream flow conditions are specified. The velocity is as- 
sumed to be constant at stations ah and de. Hence, the 
stream function varies linearly from a to h and from d to e. 
Since the angular distance from a to h and from d to e is one 
pitch angle and the flow is assumed to vary periodically 
about the axis with a period of one pitch angle, the specifica- 
tion of conditions at points a, d, and h fixes the solution. 
The value of f is constant along both the driving and trailing 
faces, with the difference between the function values being 
equal to ^ a - The choice of 0 and ±1 as the values of 
at points a, h, d, and e and along be and gf represents no 
restriction since these three basic solutions remain solutions 
to the homogeneous equation when changed by a multipli- 
cative or additive constant. 

Coefficients of i/'o, 4>i, and \p 3 in linear combinations. — 
The final solution 4' for any weight flow or rotational speed 
will be obtained from an equation of the form 

■*=Aoip a +A 1 ip 1 +A 2 t 2 +A 3 $ s (3) 

The coefficients A 0 , A,, A 2 , and A 3 are determined by the 
specification of four independent physical conditions: (1) 
the rotational speed, (2) the weight flow, (3) the location of 
the rear stagnation point, and (4) the irrotationality of the 
absolute flow. 

The coefficient Ac is determined by the rotational speed 
and is given by 

Ao=^ (4) 

Wo 

That is, Ac is the ratio of the rotational speed oj for the 
desired solution to that used in obtaining the basic solution 
0)0- 

The change in 4' across one blade passage is equal to the 
weight flow through a single passage. Therefore, 

Ai~\- AiA A 3 =M (5) 

where M is the weight flow through a single passage. 

The rear stagnation point is assumed to be at the blade 
tip. Thus, at the tip 

w s ,t — 0 
or 



This derivative is expressed in finite-difference form for the 
grid point at the blade tip and with equation (3) yields a 
linear relation in Ac, A lt A 3> and A 3 . 

The absolute flow is irrotational; therefore, if r, and r 2 
are radial stations upstream of the blade row, the following 
equation holds: 

J » 2 t r»r 

(w9.2-f-ar 2 )r 3 c#— (wb lAccr^ dd=0 (7) 
0 Jo 


where the subscripts 1 and 2 indicate values along tho lines 
r=r, and r=r 2 , respectively. If r l is chosen equal to tho 
value of r at the upstream boundary, equation (7) becomes 

X lr 

r i We, i dB=—2oyjni ( 8 ) 

since w t>1 is equal to — ur t . When the stream-function 
definition (eq. (la)) is introduced, equation (8) becomes 

X 2r d'P’ 

d6=2rwpb 2 r 3 (9) 

Equation (9) can be integrated numerically to yield a linear 
relation in Ac, A u A 2 , and A 3 . 

Equations (4) to (6) and (9) form a system of four inde- 
pendent simultaneous linear equations in four unknowns, 
Ac, A lt A 3 , and A 3 . 

Numerical method of obtaining basic solutions. — Tho 
region of solution was covered with a network of grid lines 
whose intersections form grid or nodal points, as shown in 
figure 1. 

The solution for a given set of boundary conditions of tho 
differential equation was obtained at each of these grid points 
by solving the set of linear simultaneous equations obtained 
when the differential equation is written in finite-difference 
form for each grid point. A five-point system was used in 
the finite-difference approximation of the derivatives. This 
procedure is equivalent to approximating the stream func- 
tion by a fourth-degree polynomial in the neighborhood of 
the grid point. The solution of the set of n linear simul- 
taneous equations was obtained on high-speed digital com- 
puters by the matrix method outlined in reference 12. 
Since there were four basic solutions, four sets of n simul- 
taneous linear equations were solved by this process. 

METHOD OF OBTAINING BEFINED SOLUTION IN LEADING-EDGE REGION 

Relaxation method. — Tho solution of equation (2) was 
obtained by means of relaxation techniques (ref. 13) in tho 
region outlined by the bold lines in figure 1. The residuals 
of the relaxation process were reduced to a value indicating 
unit change in the fifth decimal place of the stream function. 
This degree of accuracy is consistent with tho accuracy of 
the cross plots used to obtain the boimdary values. 

Grid. — The grid used in the relaxation process is shown in 
figure 3. The density of the grid lines is greatest in the region 
of the leading edge and decreases gradually in both direc- 
tions. This type of grid assured greatest accuracy in tho 
region of main interest. The ratio of tho grid spacing in 
either direction at a point is never greater than 2.5 and usu- 
ally less than or equal to 2.0. The intersections of tho grid 
lines are called grid points. The grid points in the interior 
of the boundaries are the points at which the numerical solu- 
tion is obtained. There were 138 interior grid points for 
this problem. 

Finite-difference approximation. — The derivatives of equa- 
tion (2) were written in finite-difference form with a three- 
point system used. In the original solution a five-point 
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system was used. Because of the fineness of the grid, a 
three-point system was considered adequate for this prob- 
lem. 

NUMERICAL EXAMPLE 

The previously outlined methods were applied in order to 
analyze the flow in a 48-inch-diameter radial-inlet centrifugal 
impeller. A description of the geometry of the impeller and 
the operating conditions for which the analysis was carried 
out follow. 

Geometry of impeller. — The impeller investigated was a 
48-inch-tip-diameter radial-inlet centrifugal impeller having 
18 blades, similar to the one discussed in references 7 and 8. 
The sharp leading edge and blunt trailing edge were rounded 
as shown in figure 1 because of practical computing considera- 
tions. The blade coordinates are given in table I. The 
solution was obtained on the surface generated by rotating 
the mean blade-height line about the axis of rotation. This 
line was approximated by the following function : 


-0.041456 
2 (r— 0.40828) 


constant 


( 10 ) 


The streamline spacing in the axial-radial plane is not known. 
Therefore, the stream-sheet thickness b in the z-direction 
was approximated by the blade height in the z-direction. 
This parameter was approximated by the following function : 

5=0.07208+1.01517e- 1 -« Mlr (11) 

The parameter X is equal to dr/dz of the stream-surface trace 
in the axial-radial plane and from equation (10) is given by 

1 0.041456 , 

X (r-0.40828) 2 { } 


Operating conditions. — Four solutions were obtained cor- 
responding to four weight flows at a tip speed of 700 feet per 
second. These four weight flows, which are the same as 
those of reference 7, are shown in the following table : 


Case 

Weight flow, 
M, lb/sec 

A 

14. 00 

B 

26. 26 

C 

32. 10 

D 

44. 00 


RESULTS AND DISCUSSION 

The results of the solutions obtained by the application of 
the previously outlined methods are presented in figures 4 
to 11, which show streamlines and constant relative velocity 
contours. The over-all solution for the entire blade passage 
as obtained by the matrix method is shown in the (a) part 
of each figure and the refined solution for the leading-edge 
region in the (b) part. Figures 4 to 11 are projections on 
the r0-plane; that is, the curvature of the stream surface in 
the axial-radial plane is neglected. The comparison with 
the results of the approximate method described in reference 
6 is made in figure 12. 


STREAMLINES 

The distribution of stream function is shown by means of 
contours of constant stream-function ratio ty/M in figures 4 
to 7 for the four weight flows investigated. The impeller 
tip speed was 700 feet per second for all four cases. 

Case A. — In figure 4 the streamlines for a weight flow of 14 
pounds per second are shown. This condition is the incipient 
surge weight flow for the experimental case (ref. 7). A large 
eddy attached to the driving face of the blade extends from 
r~1.31 to r~1.84 and almost one-third the distance across 
the passage between blades at its widest point. The major 
part of the flow is concentrated in the region near the trailing 
face, while the eddy and other relatively low-momentum 
fluid occupy half the channel. 

The inlet stagnation point occurs on the driving face of the 
blade at r~1.05, and the local angle of attack is 80°. The 
local angle of attack is defined as the angle between the 
tangents to the blade mean line and the stagnation stream- 
line. The sign convention for angle of attack, shown in 
figure 3, is chosen so that flow directed at the driving face 
results in a positive angle of attack. The slip factor, the 
ratio of the mass-averaged absolute tangential velocity of 
the fluid at the tip to the absolute tip speed, is equal to 0.874. 

Case B. — The streamline pattern for a weight flow of 26.25 
pounds per second is shown in figure 5. This weight flow 
is sufficiently high to eliminate the eddy on the driving face of 
the blade. However, a fairly large concentration of low- 
momentum air is still present, so that halfway through the 
impeller 50 percent of the fluid occupies more than two-thirds 
the available flow area. The slip factor at the impeller 
tip is 0.873, and the local angle of attack is approximately 6°. 
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(a) Entire blade passage. 
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(b) Leading-edge region. 

Figure 4. — Streamlines for weight flow of 14 pounds per second 
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(b) Leading-edge region. 

Figure 5. — Streamlines for weight flow of 26.25 pounds per socond 

(case B). 





























[iivi 


i=i;wre):~ 




Siwy 




h- — 

ID O 'f 

o m r- 

in ro w 

rf ’T. *r. 


'r ro_ o 
'f®(00 
ro f-. — 05 
id m in rt- 
f'} r<irOiO 


Biliniiifii 



1.02308 


1.01146 


39984 


8 iJ o 5 

5n O g £ 

s ? ^ 


■- l.03472r 
*o r 

§ l.02308|- 
t .0 1 1 46 j 
.999841 


mm 

"Samgg wSBS I 

■Kara 

BEBSH 

M »taa s HI S 


37656 


97656 


95332 


pounds p 


LiT:»inmr^i 


id 


Figure 7. — 
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(b) Leading-edge region. 

Figure 10. — Contours of constant relative velocity ratio TV for weight 
flow of 32.10 pounds per second (case C). 


(b) Leading-edge region. 

Figure 11. — Contours of constant relative velocity ratio W for weight 
flow of 44 pounds per second (case D). 
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Case C. — The streamline pattern for a weight flow of 32.10 
pounds per second (fig. 6) is s imil ar to that for 26.25 pounds 
per second because of the comparatively small change in 
weight flow. However, the local angle of attack changes 
considerably, from 6° for case B to —7° for case C. This 
sensitivity of local angle of attack to weight flow changes is 
quite marked over the entire weight flow range investigated. 
The weight flow for zero local angle of attack is about 29 
pounds per second. 

The flow continues to shift toward the driving face, and the 
slip factor at the tip is 0.871. 

Case D. — In figure 7, the streamline pattern is shown for a 
weight flow of 44 pounds per second. In the investigation 
reported in reference 7, this represented the maximum weight 
flow. The flow is distributed across the passage more nearly 
uniformly than in the other examples. The flow ceases to 
bo perfectly guided at i — ■■1.5, as occurred for all other weight 
flows. The slip factor at the tip for this condition is 0.859, 
so that the total variation in slip factor from surge to maxi- 
mum weight flow as determined in reference 7 is 0.015. The 
inlet stagnation point occurs on the trailing face at r~ 1.028. 
Thus, the stagnation point shifts about 0.022 foot as the 
weight flow increases from 14 to 44 pounds per second. The 
local angle of attack for case D is — 28° so that the total range 
of local angle of attack investigated was from 80° to— 28°. 

RELATIVE VELOCITY 

Contours of constant relative velocity ratio W (relative 
velocity divided by tip speed) are shown in figures 8 to 11. 
In the figures showing the entire blade. passage ((a) parts), 
the velocities near the leading edge are not shown. Reference 
should be made to the figures showing the leading-edge 
region only ((b) parts). 

Case A. — For a weight flow of 14 pounds per second (fig. 8) 
a stagnation point occurs where the eddy begins to form at 
r~1.31 on the driving face. Velocities are low along the 
entire driving face with slightly negative velocities in the 
eddy region. A rapid acceleration followed by a less rapid 
deceleration occurs on the leading edge and the trailing face 
because of the large positive local angle of attack. 

Case B. — In figure 9 the velocity contours for a weight flow 
of 26.25 pounds per second are shown. Downstream of the 
leading-edge region the velocity along the trailing face is 
nearly constant (except for a small acceleration and de- 
celeration at r~1.3 to r~1.7). In the leading-edge region, 
small local decelerations occur on both the driving and 
trailing faces with the one on the trailing face being the larger. 
Flow conditions in the leading-edge region are better for this 
weight flow than for any of the other weight flows investi- 
gated. 

Case C. — The velocity contours for a weight flow of 32.10 
pounds per second are presented in figure 10. At r~1.3 on 
the trailing face, a slight acceleration occurs followed by a 
rapid deceleration. This velocity peak is caused by the 
beginning of more rapid blade curvature at that point. 

For case C a larger deceleration occurs on the driving face 
than for case B. Decelerations are probably more serious on 
the driving face than on the trailing face because the low- 
momentum air caused by the deceleration aggravates the 
secondary-flow conditions. These secondary flows transport 
the low-momentum fluid on the driving face to the trailing 


face. This type of motion is discussed in more detail in 
reference 14. 

It was previously noted that the weight flow for zero 
local angle of attack is about 29 pounds per second. How- 
ever, zero local angle of attack is not necessarily desirablo 
for rounded leading edges such as the one considered heroin. 
In view of the velocity contours of the 26.25- and 32.10- 
pound-per-second cases, a slight positive local angle of 
attack seems desirable. For a weight flow of 29 pounds per 
second, a greater deceleration would probably occur on tlio 
driving face than that which occurs for a weight flow of 26.26 
pounds per second. This appears to be the case since the 
leading edge is shaped so that at zero local angle of attack 
the flow on the blade surfaces would be roughly symmetrical 
about the blade mean line. Consequently, it seems ad- 
visable to design the leading-edge regions in tbo manner 
suggested in reference 15 and further discussed in reference 
16. These leading-edge contours are characterized by very 
little curvature of the driving face so that flow alined with the 
driving face would produce little or no local deceleration. 

Case D. — The velocity contours for a weight flow of 44 
pounds per second are plotted in figure 11. The How along 
the driving and trailing faces is similar to that for caso C, 
because in both cases the inlet flow is directed toward the 
trailing face. The local acceleration and deceleration on the 
trailing face at r~ 1 .3 is more pronounced for this case. 

A very large deceleration occurs on the driving face just 
downstream of the leading edge. This deceleration is about 
the same size as that which occurred on the trailing face in 
case A. The deceleration is probably more serious on the 
driving face because of its contribution to the buildup of 
secondary flows. Also, at sufficiently high weight flows the 
separation following a rapid deceleration will induce choiring 
before the theoretical maximum weight flow is attained. 

MEAN ANGLE OF ATTACK 

The approximate mean angle of attack, that is, the angle 
between the mean flow direction at the inlet and the tangent 
to the blade mean line, was computed from the rotational 
speed and the average inlet velocity. The average inlet 
velocity was computed from the weight flow and the annular 
area. Two values were used for the annular area: (1) the 
total annular area with no blade blockage assumed, and 
(2) the total annular area minus the blockage caused by the 
blades. The thickness of the blades used in the latter 
computation was that at the 1.04-foot radius, which was 
approximately the radius at which maximum blade thickness 
in the tangential direction occurred. The moan anglo of 
attack across the passage at the 1.04-foot radius was also 
computed from the exact solution. These are compared in 
the following table: 


Weight 
flow, M, 
lb/sec 

Mean angle of attack, deg, based on — 

Unblocked 

annulus 

Blocked 

annulus 

Exact 

solution 

14. 00 

12. 3 

9. 2 

0 

26. 25 

. 1 

-4 6 

-6. 3 

32. 10 

-4 9 

-10. 0 

-9. 0 

44 00 

-13. 6 

-19. 0 

-16. 4 
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The sign convention for the angle of attack is such, that a 
positive angle of attack indicates flow directed toward the 
driving face of the blade. From the comparison of these 
angles of attack, it is apparent that the mean angle of attack 
is best predicted by basing the calculations on the annular 
area with blade blockage considered. The poor agreement 
between the mean angle of attack of the exact solution and 
that based on the blocked-inlet annular area at the lowest 
weight flow is probably caused by the eddy. It appears that 


inlet flow alined with the driving face results in good flow 
conditions in the leading-edge region. The blade angle of 
the driving face just downstream of the rounded leading 
edge is 57°, whereas the angle between the mean line and the 
radial direction is 62°. Thus, for case B the average inlet 
flow angle would approximately equal the driving-face blade 
angle. Flow conditions in the leading-edge region for case B 
seemed to be the best of the conditions investigated. 






(o) Weight flow, 32.10 pounds per second (case C). (d) Weight flow, 44 pounds per second (case D). 

FiaoKE 12. — Comparison of blade-surface velocity ratios obtained by approximate and matrix methods. 
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COMPARISON WITH APPROXIMATE METHOD 

Tlie blade surface velocities for the four weight flows were 
computed by the rapid approximate method of reference 6 
in order to determine the accuracy of the approximate 
method, especially in the inlet region, by comparison with 
the matrix solution. In figure 12, the surface velocity ratios 
as computed by the approximate method are compared 
with the surface velocity ratios as determined by the matrix 
method. The velocity ratios are plotted against distance s 
along the blade mean line. 

In the approximate method the “average” values of W 
and fi are values which satisfy the one- dim ensional continuity 
equation. The “average” velocity ratio is assumed equal to 
the average of the blade surface velocity ratios W ac . The 
“average” flow angle is assumed equal to p ae , the average of 
the blade surface angles, for values of 8< 1.34. For s>1.34, 



Distance along blade mean line, s t ft 


(a) Weight flow, 14 pounds por second (case A). 
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Distance along blade mean line, s, ft 
(c) Weight flow, 32.10 pounds per second (case C). 


the “average” flow angle is approximated by a parabolic 
variation (with r) between the average of the blade surface 
angles at s=1.34 and the flow angle at the tip as determined 
by the slip factor. The slope dp at /dr at s=l .34 is used as the 
third condition to determine the parabola. The approximate 
flow angle as given by the parabolic variation is also denoted 
by jS ac . It is difficult to evaluate these assumptions sepa- 
rately. In order to get some idea of their validity, the mass- 
averaged flow angle /9 a is compared with in figure 13, and 
the mass-averaged velocity ratio W m is compared with W a , 
in figure 14. 

It was impossible to obtain meaningful results in the region 
of the rounded leading edge. This failure of the approximate 
method is caused by the invalidity of the assumption that 
is equal to /3„. For values of s such that 0.07 <s< 0.6, 
the surface velocities are predicted adequately for cases B 
and C but poorly for cases A and D. For case A the failure 




Figure 14. — Comparison of mass-averaged relative velocity ratio W m with average of blade-surface velocity ratios TF a> . 
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is probably caused by the change in. the relative size of P a c 
and /3 m . At s=0.6, /3 0P is greater than (3 m but then becomes 
less than /3„ as s decreases. This trend is not followed in the 
other cases. For case D the agreement of W ac with W m 
and of ft,, with /9 ra is poorer for s<0.6 than for the other 
cases. The trends for case D, however, were the same as for 
cases JB and C. 

In the region 0.6<s<1.8, the agreement between the 
approximate solution and the exact solution is adequate 
(except in the neighborhood of s=0.8 for cases B, C, and D) 
for approximately predicting both the velocity and the veloc- 
ity gradients for all four weight flows. Such good agreement 
is surprising for case A, because the agreement between W m 
and W a , and between /S m and ft,, is not nearly so good for 
case A as it is for the other cases; the eddy attached to the 
driving face for case A is the cause of this poor agreement. 

For s)>1.8, the agreement between the approximate 
method and the exact method is poor, and ft,, is less than 
fS m for all four cases except near the tip in case A. This reverse 
in trend, together with the fact that the blade surface angles 
differ by a large amount, probably is the reason for the 
failure of the approximate method near the tip. In the 
region s>1.34, ft,, is determined by a parabolic approxima- 
tion to () m . Since ft,, and /S m agree fairly well at s=1.34 and 
at the tip, the failure of the parabolic approximation must 
be attributed to incorrect values of df}„/dr at s=1.34 or to 
the inadequacy of a parabolic approximation. 

SUMMARY OF RESULTS 

A method for the solution of the incompressible nonviscous 
flow through a centrifugal impeller (including the inlet 
region) is presented and is applied to a 48-incli-diameter 
centrifugal impeller. Solutions for the entire blade passage 
were obtained for four weight flows ranging from incipient 
surge to maximum as determined by actual impeller tests. 
In addition, these solutions were refined in the leading-edge 


region. The blade surface velocities obtained by the matrix 
solutions for the entire blade passage were compared with 
those obtained by a rapid approximate method of analysis. 
The following results were noted: 

1. A large eddy formed on the driving face of the blade at 
the incipient surge weight flow but was not present for the 
three higher weight flows. 

2. The slip factor varied from 0.874 to 0.859 as the weight 
flow increased. 

3. For weight flows of 26.25, 32.10, and 44 pounds per 
second, a local acceleration followed by a rapid deceleration 
occurred on the trailing face of the blades at a radius of 
about 1 .3 feet, that is, where the blade begins to curve more 
rapidly. 

4. The stagnation point shifted from the driving to the 
trailing face of the blade as the weight flow increased, while 
the local angle of attack varied from 80° to —28°. 

5. A small positive local angle of attack seems desirable for 
rounded leading-edge blades. 

6. Of the four weight flows investigated, flow conditions 
around the blade nose were the best for the weight flow 
corresponding to an angle of attack of —4.6° computed from 
blade speed and an upstream radial-axial velocity for which 
blade blockage had been taken into account. 

7. The rapid approximate method was adequate in the 
inlet region just downstream of the leading edge for weight 
flows of 26.25 and 32.10 pounds per second. These weight 
flows corresponded to mean angles of attack of— 5.3° and 
—9.0°, respectively. 

8. The mean angle of attack was best predicted by basing 
the approximate computation on the weight flow, the tip 
speed, and the annular area minus the blockage of the blades. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 4 , 1955 
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APPENDIX 


SYMBOLS 


The following symbols are used in this report : 

An, A u coefficients in eq. (3) 

Aj, A3 

b stream-sheet thickness in 2 -direction, ft 

M weight flow through single passage, Ib/sec 

r radial distance, ft 

s distance along blade mean line, ft 

ll r ratio of relative velocity to tip speed 

v) relative velocity, ft/sec 

2 axial distance, ft 

/S flow angle, deg 

6 angular coordinate in relative system, radians 

X slope of trace of stream surface in axial-radial 

plane 


p 

'S' 

Mu 

0) 


fluid density, Ib/cu ft 
stream function, eqs. (1) 
basic solutions 

angular velocity of impeller, radians/sec 


Subscripts : 
av 


m 

r 

t 

6 

0 

1,2 


approximation to averaged value used in approx- 
imate method 
mass-averaged value 
component in radial direction 
impeller tip 

comporent in tangential direction 

value in basic solution of eq. (2) 

conditions along r=rj and r— r*, respectively 
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TABLE I.— MODIFIED BLADE COORDINATES 


Driving face 

Trailing face 

r,ft 

9, radians 

r, ft 

9, radians 

1. 0405 

0. 34256 

1. 0285 

0. 34890 

1. 0521 

. 34890 

1. 0352 

. 37797 

L 0696 

. 37224 

1. 0405 

. 38604 

L 0740 

. 37797 

1. 0518 

. 40705 

L 1190 

. 43612 

1. 0696 

. 43273 

L 1278 

.44633 

1. 0971 

. 46520 

1. 1711 

. 49427 

1. 1278 

. 49939 

L 2324 

. 55242 

1. 1611 

. 52335 

L 2441 

. 56242 

L 2144 

. 58150 

L 3095 

. 61057 

1. 2441 

. 60413 

L 3604 

. 63996 

1. 2974 

. 63965 

L 4262 

. 66872 

1. 3526 

. 66872 

L 4767 

. 68476 

L 3604 

. 67227 

L 5930 

. 70726 

1.4283 

. 69780 

1. 7093 

. 71497 

1. 4767 

. 71113 

1. 8256 

. 71601 

L 5600 

. 72687 

L 9419 

. 71601 

1. 6930 

. 73041 

2. 0000 

. 72080 

1. 7093 

. 73583 

Z 0123 

. 72687 

1. 8256 

. 73416 



1. 9419 

. 73312 



Z 0000 

. 73145 



Z 0123 

. 72687 




